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Abstract 

We study the essential spectra of formally self-adjoint elliptic systems on doubly periodic 
planar domains perturbed by a semi-infinite periodic row of foreign inclusions. We show 
that the essential spectrum of the problem consists of the essential spectrum of the purely 
periodic problem and another component, which is the union of the discrete spectra of model 
problems in the infinite perturbation strip; these model problems arise by an application of 
the partial Floquet-Bloch-Gelfand transform. 

Keywords: Spectral bands, double periodic medium, unbounded periodic perturbation, 
open waveguide 
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1 Introduction. 

1.1 Preamble. 

Composite materials are used extensively in the modern engineering practice. They are math¬ 
ematically interpreted, modelled and treated by applying the Floquet-Bloch-Gelfand-theory 
(FBG-theory in the sequel) for elliptic spectral boundary value problems with periodic coeffi¬ 
cients in periodic domains. This approach has led to many important theoretical results and 
applications in topics like homogenization, diffraction in waveguides, band-gap engineering etc. 
The usnal setting for the theory concerns purely periodic media (periodic coefficients, periodic 
domains), though certain types of perturbations may be allowed as well. In this paper we 
introduce and examine quite novel type of perturbations of a double-periodic medium with 
semi-infinite rows of foreign inclusions as depicted in Fig.dl a), b). The perturbation may in¬ 
fluence the essential spectrum of the problem, and the description of this spectrum becomes 
the main goal of onr paper. Such perturbations of periodic lattices have not been the sub¬ 
ject of thorough mathematical investigation, but in the physical literature they are related, for 
example, to defected photonic crystals, cf. [m Gh.7]. 
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Figure 1: a) Domain D (grey) including the open l-shaped waveguide (dark grey), b) Strip 
(dark grey), which coincides with the open waveguide in the special case = 0. 
Dotted lines indicate in a) the cell and in b) the cell . 
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Figure 2: a) Purely periodic plane b) Plane with infinite perforated strip (dark grey) 
and the strip 11 (between dotted lines). 


In the case of homogeneous media, foreign inclusions like inhnite or semi-inhnite strips are 
called open waveguides; we accept the same terminology in the periodic case. The physical 
meaning of the notion is clear, see mu for acoustics and [HI III1EI1E3 for similar defects 
of periodic media in solid state physics and optical systems. Our approach can readily be 
generalized for different shapes of insertions, see Fig.[3l These types of perturbed periodic media 
may also appear in applications like composite materials, because of improper manufacturing or 
also a specific feature created on purpose. For technical simplicity we will deal with particular 
open waveguides as depicted in Fig.[Tl but in Section |4] we describe minor modifications of 
our approach for shapes like in Fig. [3] as well as other generalizations, including smoothness of 
coefficients and the boundary. 

We will consider a general elliptic spectral problem, see (I1.18h - (I1.19D . in the perturbed 
domain D C Fig.[TJ The main result of the paper. Theorem [7] in Section 13.51 contains 
the following statement: the essential spectrum aes (-4) of the main problem is the union of 
two subsets, the first of which is the essential spectrum aes (-4'^) of the problem on the purely 
periodic domain see Fig. [21 a), and the second one, denoted by which is caused by the 
perturbation of the domain, the open waveguide see Fig.[Tl The subset a'^ equals the union 
of the discrete spectra of a family of model problems on a domain Fig. [21 b), which is related 
to We will discuss the structure of aes (*4) hr more detail in Section [LTl after presenting 
the notation and dehnitions. 
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1.2 Purely periodic medium. 

We consider the geometry of the unperturbed perforated plane and the related spectral bound¬ 
ary value problem. Let Q = Qi C be the unit open square, Qa = (0, a) x (0, a) , and let 
be an open set in the plane with closure \J C Q and a smooth boundary 

d'd of Holder class where <5 > 0. Here is a not necessarily connected set describing 

perforation; neither the case = 0 is excluded. We define the periodicity cells (outlined by 
dotted line in Fig.[Il b) 

uj^ (a) = jx = {xi,X 2 ) : (xi - 01 , 3:2 - 02 ) G := Q \ (1.1) 

where a = (oi, 02 ) G is a multi-index and Z = {0, ±1, ±2,...} . The perforated plane is 
covered by the periodicity cells (|l.ip . so it satisfies 

\ IJ (1,2) 


where (o) denotes a translation of similarly to (jl.ip . Assuming uj^ to be a domain, we 
observe that includes the boundary strip 

{x G Q : dist(x,(9Q) < d} (1.3) 

with some d > 0 and, thus, is a domain and in particular a connected set as well. 

We consider the spectral problem 

(x, V) (x) = Au° (x), X G (1.4) 

N°(x,V)u°(x) =0, xg9H°, (1.5) 


and its variational form 

(AOzi(V)uO,Il(V)uO)^o = A(uO,uO)^o, e . (1.6) 

Concerning the notation, by A we understand a spectral parameter and u = (ui,..., tt^)^ is a 
vector function, realized as a column, so that T stands for transposition. The matrix differential 
operators and in (m and (II.5p take the form 

(x, V) = D (-V)^A° (x) D (V), (1.7) 

iV° (x, V) = D{v (x))^A° (x) D (V), (1.8) 

while V = grad, v = (z^^,z/^)^ is the unit outward normal vector on and is a matrix 
function of size m x m with C^A.gjxLooth entries; AP ( x ) is Hermitian, positive definite and 
1-periodic in x G (it is convenient to define the coefficient matrix in the intact plane). 
Furthermore, D (x) is a matrix function of size mxn and linear in x = (xi, X 2 ), and D(0, 0) = 
Onxm is the null matrix. The substitution xj djdxj gives an mxn-matrix D (V) of first order 
differential operators with constant complex coefficients, while L^(x, V) (respectively, N^{x, V)) 
is a formally self-adjoint differential matrix operator in divergence form (resp. the corresponding 
Neumann condition operator); the bar in (11.71) and (II.8p indicates complex conjugation. Finally, 
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( , )qo is the natural scalar product in the Lebesgue space (^*^)) and is the Sobolev 

space with standard norm 

\\w,H^ (0°)|| = (\\Vw,L^ (0°)||^ + ||u;;L2 

The last superscript n in the integral identity (II.6p shows the number of components in the test 
function v = (ui, ...,Vn)'^ , however, this index is omitted in the notation of norms and scalar 
products. In this way, the left- and right-hand sides of (jl.hp involve scalar products in 
and respectively. 

We assume that D (x) is algebraically complete [40]: there exists a € N = {1,2,3,...} 
such that for any row p = {pi,...,pn) of homogeneous polynomials in ^ = (^ 1 ,^ 2 ) € 1^^ of 
common degree q'> qd, one can find a row of polynomials q = {qi, ...,qn) satisfying 

P (0 = 9 (0^(0, vcgm2. (1.9) 

According to |40l §3.7.4], this assumption yields the Korn inequality 

||u; (0°) f < CD,^o (||D (V) u; (a;°) |p + ||u; (a;°) f) (1.10) 

and, hence, the sesquilinear Hermitian positive form on the left in (|1.6I) is closed in 

(n'^)"’. Moreover, the operator is elliptic and the Neumann boundary operator covers 
it in the Shapiro-Lopatinskii sense everywhere on (see, e.g., |35l Thm. 1.9]). 

Owing to the above-mentioned properties of a®, the problem (II.4p . (jl.5p is associated with a 
positive self-adjoint operator in (O*^)”' with the differential expression {x,'V) and the 
domain 

V (M°) = {u° G (0°) : iV° (x, V) (x) = 0, x G 511°} ; (1.11) 

see lU Ch. 10] and [44] Ch. 13]. The description of the spectrum a (-4.°) is well known and will 
be presented in Section [2j 

We emphasize that the results of the paper remain valid for other types of boundary con¬ 
ditions, in particular, for the Dirichlet conditions, cf. Section im iil.lznl. A description of 
all admissible boundary conditions can be found in m and [Ml §!]■ Moreover, the 
smoothness of the boundary was assumed in order to have the elementary formula (jl.lip for 
the operator domain V (M°) and to simplify the technical computations. The results of our 
paper hold true in the case of uniformly Lipschitz boundaries, but the proofs would require 
small modifications, cf. Section 

1.3 Periodic medium with open semi-infinite periodic wavegnide. 

In this section we describe the geometry of the open waveguide and the fnll spectral problem. 
Let •Cc7+ be the semi-strip {x : xi > 0, |x 2 | < h} with some h G N (overshaded in Fig.[Tlb). In 
the rectangle = (0,1) x {—h, h ), we introduce an open set with a smooth boundary d'd^ 
and closure C Q^. We define a semi-infinite row of holes or inclusions as depicted in Fig.[Tl a) 
or b), respectively: 

uj^ (ai) = {x : (xi — «!, X 2 ) G , ai G Nq = NU {0} , (1-12) 

\ 51 , (ai) = {x : (xi - oi, X 2 ) G 5^} . 
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The cell is outlined in Fig. [Ha), by dotted line. We also introduce a smooth Hermitian 
matrix function in entries of which are supported in vj^ and become 1 -periodic in xi 
inside the semi-strip {x € w~^ : Xi > i?} , i? € N, 

{x) = {x) for xi > i?, (xi ± 1 , X 2 ) = (xi, X 2 ). (1-13) 


Furthermore, 

= (J #(q;i), = IJ uj^{ai), 

ai€No oiSNo 


and the sum 

A (x) = A^ (x) + A^~^ (x) 
is assumed to be positive definite in O, where 


(1.14) 


(1.15) 


fl = (fl°\ti7+) (1.16) 

In other words, we make a perturbation of coefficients and boundary inside the semi-strip w~^, 
see Fig.[TJ For instance, one may suppose that A^^ is the null matrix Omxm and 7 ^ 0 , 
= 0, h = 1, which means filling in all holes inside ro"*", cf. Fig.[Tlb). Vice versa, in the case 
= 0 , 7 ^ 0 , one perforates the plane with a semi-inhnite row of holes, see Fig.dja). 

Even in the case of absence of holes we still call (|1.14l) the perforated strip, and the perforated 
plane in (IL21) can also contain no holes. 

Notice that the Korn ineqnality 


u;H^{n)\\'^ <cd(w°,w^) (\\D {V) u] L'^ {n)\f + \\u] L'^ (1.17) 


is still valid and can be derived by snmming np ineqnalities of type (ll.lOp in the cells 0 ;*^ (a) 
with a («! > 0, —h < a 2 <h excluded) and (ai) with ai G No- 

Replacing A^ (x) with (|1.15p in (|1.7I) still gives an elliptic and formally self-adjoint matrix 
operator L (x, V^,). The same change in (|1.8I) yields the Neumann bonndary condition operator 
N (x, Vj;), where n is regarded as the outward unit normal on 911. In the domain (jl.I 6 p we 
consider the spectral problem 


L (x, V) u (x) = Au (x), X G n, (1-18) 

V (x, V) n (x) = 0, X G 9n, (1-19) 

and the corresponding integral identity 

a{u,v) := {AD{V)u,D{V)v)^ = X{u,v)^, Vu G 77^ (Sd)” . (1.20) 

Since a remains as a closed positive Hermitian form in (f7)” , the variational formulation 
(jl. 20 p of the problem (ll.lSp . (ll.lOp supplies it with a positive self-adjoint operator A in (fl)"’ 
with the differential expression L (x, V) and the domain 

V (^) = {u£H^ (0°)” : N (x, V) n (x) = 0, x G 911} ; (1.21) 

see again [U Ch. 10] and |44[ Ch. 13]. This notation is quite similar to the one nsed in Section 

[L21 
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The main goal of the paper is to describe the essential component Uf-g (^) in the spectrum 


of A, 

a (Al) = adi (Al) U aes (A). 

(1.22) 

We emphasize that in general 

aes (Al) A (^es (Al°) , 

(1.23) 

and moreover, the discrete spectrum of A.'^ is empty, thus, 



O 

b 

II 

O 

b 

(1.24) 

We will identify the difference 

CTad (Al) = aes (Al) \ a (Al°) , 

(1.25) 


but leave aside two interesting and important questions. First, we are not able to describe 
completely the component (.A) in a general perturbed problem (I1.18p . (I1.19p . although, of 
course concrete examples of isolated and embedded eigenvalues in a^o (-A) can be constructed in 
scalar problems. Second, the existence or absence of the point spectrum fjpo (eigenvalues 
of infinite multiplicity) in the purely periodic problem (II.4p . ()1.5I) remains unknown; note that 
this question is answered in the literature only for particular scalar problems (see, e.g., papers 
[l6l[T8l[2l] and books [l3l[28]). Notice that a^o can be included in a^o (*A), but the latter 
stays unknown, too. 

1.4 Discussion on the main result 

In the case of the purely periodic plane Fig.[2la), the spectrum Oes of the problem 
(|1.4p - (|1.5h has representation as a union of spectral bands (see (I2.15p . ()2.16p . below), which 
is a well-known consequence of the FBG-theory; we refer here to [271 [281 US]. Consider for 
a moment the domain of Fig.O b) with foreign inclusions or holes, which form a periodic 
row, inhnite in both directions. Then, the spectrum may be different from the purely periodic 
case, and we denote by the essential spectrum of the problem on 12^ (cf. (|3.4p . below). 
Analogously to [ITIEIIES], the increment \ aes can be detected by performing the 
partial FBG-transform in xi-direction (Section 13.ip and investigating the kernel of the model 
problem in the perforated strip 11 (separated by dashed lines in Fig.|2l b); cf. (j3.1ip . (j3.9p ). This 
problem depends on the Floquet parameter (, and, for certain values of the spectral parameter 
A € (0,-(-oo) and ( € [0,27r), it can have a solution in the Sobolev space 27^(11)”'. Such values 
A form the point spectrum Upo (aI** (C)) of the model problem (I3.30p for the operator (C)- 
Our main result in Theorem [7] says that the following formula holds true for the problem 
(ll.lSp . (I1.19P in the periodic plane with the immersed l-shaped open waveguide, see (ll.lGp and 
Fig.dlb): 

CTad (Al) = cr** \ aes (aI°) = |J apo (a^ (C)) \ CTes (aI°) . (1.26) 

CG[ 0 , 27 r) 

The last set in this formula requires some comments. First, embedded eigenvalues in aem {A^ (C)) 
<7po {A^ (C)) \ ^di (aI^ (C)) live inside the essential spectrum aes (aI** (C))) which in turn is con¬ 
tained in aes (aI*^) (compare (|3.32l) with (I2.15p l. Second, aes {A^ (()) depends on ( and therefore 
some points of the discrete spectrum adi (aI^ (C)) may fall into aes (aI** (CO) with / (. In any 
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case none of the indicated points in Upo {A^ (C)) stays in the increment component (11.251) of 

^es (•^) • 

We emphasize that, contrary to the case of the lacking periodicity of the domain 
prevents a direct use of the partial FBG-transform, hence, the proof of (|1.26l) requires improved 
mathematical tools. The new procedure of our paper involves the construction of a parametrix 
and singular Weyl sequences in order to describe, respectively, the regularity field and the 
essential spectrum of the operator A. 

1.5 Cranked and branching open wavegnides 

Fig. [3] shows open waveguides of the shape of the letters V, X and Y, Z. They appear due to 
perforation and perturbation of coefficients in overshadowed joints of semi-strips. In Section 
m we explain how the results of Section [3] for the l-shaped case, Fig.[Tl a), b), can be readily 
adapted to these cranked and branching open waveguides. Here, skewed branches of the V- and 
Y-shaped waveguides must maintain the periodicity and thus the tangent of tilt angles has to 
be rational number. We do not know a formula for the essential spectrum in the irrational case. 

We also emphasize that for a clear reason, no relevant perturbation in a disk B/j = {x : r = \x\ < R} 
with radius R > 0 can affect the essential spectrum of the boundary value problem (|1.18l) , (11.1911 . 
Moreover, assume that A [x) is a matrix decaying together with its derivatives at infinity as 
(1 -|- , 5 > 0, and that the coefficient matrix ()1.15p 

A{x) + A{x) (1.27) 

still keeps the above mentioned basic properties of A. This replacement of the coefficient matrix 
does not change aes- 

All these generalizations and some others will be commented in Section HI We have chosen 
the very particular open waveguide in Section 11.31 in order to simplify the presentation, to 
illuminate the main points of our approach and to avoid unimportant but cumbersome technical 
details. 

1.6 Structure of the paper 

In Section [2] we recall generally known information on the purely periodic case which will be 
used later on. The main interest is focused on the model problem (j2.1ip in the periodicity cell 
w, which is obtained using the FBG-transform |19] . The open periodic semi-infinite waveguide 
will be considered in Section [3l where we apply the partial (TGil)-transform to formulate 
another model problem (j3.30p in the perforated inhnite strip H with periodicity conditions on 
its lateral sides. 

The spectra of those two model problems form the essential spectrum Ges (.4) of the original 
problem (ll.20p . To verify the corresponding formulas (|3.38l) and (|2.16l) . (I3.39P we first present 
two types of singular Weyl sequences for A and on the other hand construct a right parametrix 
for the formally self-adjoint problem (ll.lSp . ()1.19p . This is the most involved part of our paper. 

To that end, we follow [33] and also |39p §3.4], and study an operator family for a second model 
problem in the weighted Sobolev spaces (Kondratiev spaces) kF^ (H) , leading to important 
conclusions on exponential decay properties of the solutions in the strip H. Finally, we glue 
the parametrix (|3.46p from solutions of the model problems with the help of appropriate cut-off 
functions. The parametrix enables to prove that, for any A G M+ outside the union of sets ()2.16p 
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Figure 3: V,X and Y, Z-shaped joints of open waveguides (dark grey). 


and (I3.39|) . the operator of the inhomogeneous problem (jl.lSjl . (I1.19h . cf. (I3.4njl . is Fredholm in 
the Sobolev-Slobodetskii spaces; therefore such points A form the intersection of the regularity 
field of A with the semi-axis M+. This completes the proof of Theorem [71 the central assertion 
in the paper. 

We start the last section of the paper by describing several concrete problems in acoustics, 
elasticity and piezoelectricity, to which our theory may apply. However, as mentioned above, 
the original exact formulation of the spectral problem was simplihed in several aspects, so 
we comment in the next subsections on certain supplementary issues in order to obtain more 
generality for further interesting physical applications. We finish the paper with Section 14.51 
where we present small modihcations of the parametrix to be applied to semi-bounded open 
periodic waveguides in the shape of the letters V, X and Y, Z as in Fig.[3l 

2 Spectrum of the purely periodic problem 

2.1 Floquet-Bloch-Gelfand-transform 

The spectrum of the operator Xl® in the purely periodic domain can be studied with the 
help of the FBG-transform. This will lead to the formula (I2.16p . the main object of Section [2l 
The FBG-transform is defined by 

u 1 -^ u {x; rj) = -^ ^ {x + a), x € oj^, (2.1) 

aGZ2 

and it establishes the isometric isomorphism 

(0°) ^ (Y; (a;°)) (2.2) 
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(cf. [T9] and, e.g., 


, where r} ■ x = rjiXi + rj 2 X 2 , 


Y = {r]= {r]i,Tj 2 ) : G [0,27r)} , Y = Q 2 ,r, 


(2.3) 


and (Y; *B) is the Lebesgue space of abstract functions in r/ G Y with values in Banach space 
IB and the norm 


( rZTT rZTT 


1/2 


Moreover, the mapping 




(2.4) 


(2.5) 


is an isomorphism, too. Here, subspace of vector functions U G 

satisfying the periodicity conditions 


dU dU 

U {0,X2) = U {1,X2) , ^ (0,X2) = ^ (1,X2) , X2G(0, 1) 

H(xi,0) = C/(xi, 1) , ^ {xi, 0 ) = ^ {xi,l), xiG(0, 1) 

0X2 0X2 

In what follows we shorten the notation dU/dxj to dj U. 

The inverse FBG-transform is given by 

1 /•27r r2n 

U = u>-B u{x) = — / / (x — [x] ; r]) dr]idr]2, 

2vr Jo Jo 

[x] = ([xi] , [X2]), [t] = maxjr G Z : r < t} . 


( 2 . 6 ) 


(2.7) 


2.2 The model problem on the periodicity cell 

Owing to ()2.ip . we have 

{x-r]) = {di+ir])'^^ {d2 + ivT^v{x]T]), V“ = (2.8) 

and thus the FBG-transform (|2.ip converts the problem (II.4p . pi.5p into the following problem, 
depending on the parameter r/ G Y, in the periodicity cell co^ , 

{x,V + ir])U {x;r]) = A{ri)U {x;r]), x e lu^, (2.9) 

{x,V + ir])U {x;r]) = 0, x € = div^ Ci Q, (2.10) 

together with the periodicity conditions (12.6p on the exterior part dQ of the boundary of the 
cell uj^. The variational formulation of problem (12.9p . (12.101) . (12.61) amounts to finding a number 
A{rj) ^ C and a non-trivial vector function U (•; ?/) G such that 

{A^D (V + iv) U,DiV + ir^) V)^, = A (r?) (U, H)^o , VH G (a;°)" . (2.11) 

In view of the compact embedding (o;*^) C (w®) in the bounded domain cj^, the spectrum 
(7° of the variational problem (12.111) and boundary value problem (|2.9p , (|2.10l) , (12.6p is discrete 
and forms the unbounded monotone sequence 

0 < Ai [r]) < A2 (7?) < ... < Afc (7?) < ... +00, ( 2 . 12 ) 
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where multiplicities are counted. According to the general results of the perturbation theory 
for linear operator^, the functions 

Y3rj^Ak{v) (2.13) 

are continuous, see for example [231125| . Moreover, they are 27r-periodic in rji and because 
for any eigenpair {A ( 77 ), U (x; rj)} of the problem (12.9|) . (I2.10p . (12. 6 p at some 77 G Y, 

{A(77),e±2™p[/(x;77)}, p=l,2, (2.14) 

remains an eigenpair of the same problem but at 77 ± 27re(p) where e(p) = (Si^p, 62 ,p)~^ G i 
the unit vector of the 77 p-axis. 

The above mentioned properties of functions (|2.13p ensure that the spectral bands 

= {Afc ( 77 ) : 77 G Y} , /c G N, 

are bounded connected closed segments in M_|_. The formula 

a{A^)=aes (- 4 °) = 

fceN 

for the spectrum of the problem (jl.4p . (II.5p is well-known, see, e.g., [271 ill [28], but we briefly 
comment on its proof in Sections 12.3112.41 since we will need some of these arguments later. 

2.3 Unique solution of the inhomogeneous problem 

As regards to (12.161) . we now prove the inclusion a (M'^) C Let us consider the 

boundary value problem 


(2.15) 


(2.16) 


{x, V) (x) — (x) = (x), X G 

(x, V) 77 ° (x) = < 7 ° (x), X G 


(2.17) 


with the data 


f G (flO)'^ ^ ^0 g jji /2 ^2.18) 

and a fixed parameter A = such that 

A° i Bl, yk G N. (2.19) 

In (I2.18p . (cin^) stands for the Sobolev-Slobodetskii space of traces with the intrinsic norm 

{dn^) =inf{||G°;F^ ( 0 °)|| : = 5 ° on (90°} . ( 2 . 20 ) 

This norm is equivalent to the following one: 

\ 1/2 

Q / N I 2 dSxdSy \ 


f||(7°;L2(aO°)||^+ / [ |,°(x)-(7°(2/)| 

y Jdno Jdn°{x) 


\x - yr 


Here, dsx is the arc length element on 90° and 90° (x) = {y G 90° : |y — x| < /} while I G 
( 0 , -|-oo) can be fixed arbitrarily. 


quadratic pencil easily reduces to a linear non-self-adjoint spectral family. 
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Clearly, the mapping 

{ny {lV - A%0,iV%0} e {ny X ^Q^oy 


( 2 . 21 ) 


is continuous for any A^, but in the case (12.191) it becomes an isomorphism. Indeed, the FBG- 
transform ()2.ip turns ()2.17l) into the parameter-dependent problem 


L'^ {x,V + ir])u^ {x;r]) — {x;ri) = f {x;ri), x G a;°, 
(x, V + ir])u^ {x;r]) =g^ (x; ry) , x G 7*^, 

with the periodicity conditions (|2.6p . The right-hand sides meet the estimate 


( 2 . 22 ) 


/ L2 ( 0 , 27r; y))f + jj? (o, 27r; 77^2 (^ 0 )^ 

< c (^||/°; ( 0 °) II' + || 5 °; ^ 1/2 (ailO) |Q , 


(2.23) 


while the necessary information about the Sobolev-Slobodetskii space is provided by the iso¬ 
morphisms (|2.2I) . (12.4p and the formula (|2.8I) for derivatives. By the assumption (|2.19l) . the 
problem (I2.22p . (|2.6p has for any r/ G Y a unique solution denoted by 


u°(x;r/) = (A°;r/) {/°(-;r/) ,?°(-;h)} 


(2.24) 


and estimated as follows: 


y) f <c(jf\-,ri)-,L^ y) II' + 11? ° (•; r?) ; ( 7 °) ||') . (2 


.25) 


Since the constant C does not depend on ry G Y, it suffices to apply the inverse FBG-transform 
(12.7p and to derive from (I2.25P and (|2.23l) the inequality 

\y;H^yf<c(^\\f-,L^n°)\f+ {dn°) (2.26) 

for the unique solution G (ll^)"’ of the problem (|2.17l) with fixed parameter (|2.19l) . Thus, 
mapping (I2.21h with this A*^ is indeed an isomorphism, which in particular means that A'^ belongs 
to the regularity field of the operator Al® in (jl.lip . This coincides with the resolvent set of 
because the discrete spectrum adi is evidently empty. 


Remark 1 Dealing with the Sobolev-Slobodetskii norms in 112.2S\} . f2.25]) and \2.26\) . it is much 
more convenient to use the definition 112.20\) and some extensions G of g and G (•;??) ofg{-\if} 
such that 


g-,H^/^dn°) < 2 ||G; 771 (110)11, 
g- ( 0 , 27r; 77^2 (/)) jj < 2 ||G; (o, 27r; 77^2 (/)) 
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2.4 The singular Weyl sequence 

We next show that a (.4°) D Let us assume 

A° G B'I for some A: G N (2.27) 

so that there exist rp and (•; G such that {•\Tp ')} is an eigenpair of 

the problem (I2.9I1 . (|2.10l) . (12. 6 p with A [rp) = A°. By a direct calculation, one easily deduces 
that the Bloch wave 

{x) = P^°-^U^ {x]rP) (2.28) 

satisfies the differential equations (jl.4p and the boundary conditions (jl.5l) . although it does of 
course not fall into the Sobolev space (12*^)"'. However, (I2.28P is useful for constructing a 
singular sequence hr V (aI°) C for the operator at the point A*^, namely 

a sequence with the following properties: 

1 ° \\u^^;L^ {n^)\\ = l; 

2° —r 0 weakly in L? (l2‘^)”' as j +oo; 

3° (H°) II —)■ 0 as j —+oo. 

To dehne the entries of this sequence, we introduce the plateau function 

Xi (A) = X (A - 2^) X (2^+^ - t) , (2.29) 

where x £ (M) is a cut-off function such that 

X (t) = 1 for t > d, X (A) = 0 for t < 0, 0 < x < 1 (2.30) 

and d > 0 is taken from ()1.3p : therefore, in the vicinity of each component 7 *^ (a) = d'd^ (a) of 
the boundary the two-dimensional plateau function 

(x) = Xj (xi) Xj {X 2 ) (2.31) 


becomes a constant, either 1 or 0 . 

We set 

{x) = Xj{x)u^{x), u°^(x) = ||u°^; 7.2(12°) ||■\o^■(x). (2.32) 

The above specification of Xj shows that 

(x, V) (x) = Xj (x) {x, V) (x) = 0, x G (912°, 

and, hence, G V (aI°) . 

The property 1° clearly holds true. Furthermore, the weak convergence to 0 in 2° occurs 
at least along a subsequence of indices j G N, because, by (I2.29P and (I2.3ip . supp (u°-^) n 
supp (u°^) = 0 as j / /. It remains to verify the property 3°. Recalling (I2.32p and (|2.28l) . we 
have 

||u°^;L 2 (o0)||2 > f \U^ [x-,ip)f dx = (2.33) 

aeZ2n[2J+l,2J+l-l] 

= ( 2 ^+^ - 2 ^ - 2)^ ||[/° (•; T/°) ; (^^0^ ||2 > ^ 022 ^ 
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Figure 4: Supports of the cut-off functions. 


with Co > 0 and j > 2. At the same time, we obtain 




where stands for the commutator of the differential operator (|1.7I) and the multipli¬ 

cation operator with Xj. Owing to definition (I2.29p - (l2.3ip . the plateau function (I2.3ip varies 
only inside the union of four rectangles of size d x 2 ^ and the common area O (2^) (see the 
overshaded frame in Fig.HJ a)). Hence, due to periodicity in (|2.28p . we arrive at the inequality 


(oO)lf <ci2^ 


which together with (I2.33P and (|2.32l) prove the relation 

II ; (0°) II < C2-^/^, C = ci/cq 


as well as the property 3°. 

By the Weyl criterion (cf. [H Th. 9.12] or [321 Th.VII.12]), the point (I2.27P lives in the 
essential spectrum aes (“4*^) ■ This and the material of Section [2.31 confirm the formula (I2.16p . 

3 Spectrum of the open waveguide in periodic medium 

3.1 Partial Floquet-Bloch-Gelfand-transform 

Our aim is to apply the partial FBG-transform to detect the effect of the open waveguide to 
the essential spectrum of the problem (|1.18l) - (|1.19l) . Due to the lack of periodicity, this cannot 
be done directly in 0, hence, we introduce and study in Sections I3.1H3.3I the problem in the 
domain Dt*, Fig. [21b). The results of Section ISTHl will be applied in Sections I3.4H3.5I to the 
original problem, which leads to the proof of the main result, Theorem [71 

Similarly to (|1.14p - (jl.l6p we introduce the infinite periodically perforated strip (over¬ 
shaded in Fig. [21 b), 

= ro \ |^a;^(ai), zu = M x (—/i,/i), (3.1) 

aieZ 
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and the positive definite Hermitian matrix 

(x) = (x) + (x) (3.2) 


which happens to be 1-periodic in the variable xi. In the domain 

= (0° \zu)un^ (3.3) 

we consider the anxiliary boundary-value problem 

(x, V) (x) — X'^u^ (x) = f'^ (x), x G (3.4) 

(x, V) tt** (x) = < 7 ^ (x), X G dn'^, 


where is a fixed parameter and the operators L^, are given by (ll.7p . (|1.8I) with the change 

The domain (13.3p is also 1-periodic along the xi-axis, that is 

= |x : (xi ± 1,X2) G , (3.5) 

and we define the perforated strip bounded by dashed lines in Fig.[2lb), 


n = 


|x G n** : xi G (0,1)| . 


(3.6) 


The periodicity observed in (j3.2|) and (13.3|) allows us to apply the partial FBG-transform, 
see e.g. 


w 


^ 1 -^ [/** (x; C) = ]— e _|_ q,^. a; 2 ) 


aiSZ 


which establishes the isomorphisms 


( 0 , 27 r; (H)) , (o#) « (O, 27 r; (H)) ( 3 . 8 ) 

where the first one is isometric (cf. (j2.2l) i. while Flperf) (n)"’ is the subspace of functions U'^ G 
(n)"' satisfying the periodicity conditions on the lateral sides of the perforated strip (j3.6p 

t/“( 0 ;x 2 ) = C/“(l;x 2 ), ( 0 ; X2) = ( 1 ; X2), X2 G M. ( 3 . 9 ) 


The inverse partial FBG-transform is given by 


n“(xi,X2) = ^= / e*^*il 7 ** (xi - [xi] ,X2;C)dC, 

V 27r Jo 


(3.10) 


see [33] and cf. ()2.ip . (12.7p . 

The FBG-transform (j3.7p applies to the problem (13.4p and converts it into the parameter- 
dependent problem 

{x,di +1^,82)11'^ {x\if) — X^U'^ {x] if) = F'^ {x\if), X G n, (3.11) 

(x, di + iC, 82 ) t/** (x; rj) = G'^ (x; rf) , x G F, 


with the periodicity conditions (13.91) . Here, G M_|_ is fixed, are the FBG-images of 

and 

F = {x G an : 0 < xi < 1} (3.12) 

is the interior boundary of the perforated strip H. 
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3.2 Second model problem in weighted spaces 


To study the problem (13.111) , (I3.9|) , we introduce the weighted Sobolev space (II) (the expo¬ 
nential Kondratiev space [26]) as a completion of the linear set (ll) (infinitely differentiable 
functions with compact supports) with respect to the norm 

||u; (n)|| = {\\v\L} (n)ll' + ||vu; l} (n)|f + ||u; l} (n)|f)'^', (3.13) 


where is the family of all second-order derivatives of u, /? G M is a weight index and (II) 
stands for the weighted Lebesgue space, 


hL} (n)|| 


g/3|x2|^.^2 


(3.14) 


Notice that IT| (II) consists of all functions v G (ll) with finite norm (I3.13|) . The norm 

modified by omitting Vu; (II) on the right in (I3.13P remains equivalent to the original 

one. For /3 = 0, we have Wq (II) = (II), but in the case /3 > 0 (/3 < 0) the Kondratiev 
space includes functions with an exponential decay (growth) at infinity with decay (growth) 
rate controlled by the weight index. By IT| (II), we understand the subspace of functions 

subject to the periodicity conditions (13.9p . and (T) is the weighted Sobolev-Slobodetskii 
space with the intrinsic norm 




inf{||K;IT^i(n)|| 

(||VK;L2(n)|f + ||K;L2(n)|f)'^' : K = uonr}. 


(3.15) 


1/2 

We emphasize that W^ (T) does not require periodicity conditions because (|3.12ll includes the 
interior part of the boundary. 

Owing to definitions (I3.13l) - (l3.15p and formula (13.Sp . the partial FBG-transform establishes 
the isomorphisms 


L} (O#) (0, 27r; (B)) , W} (o#) (O, 27r; (B)) , (3.16) 

(^0, 27r; (f«)) . 

The problem operator for (13.lip . (13.9p , 

(n)"9B#^r^(A«,c) 

= {l# (x, Si + iC, ^ 2 ) - A», iV» (x, + ic, <92)} G L} (B)” x w]'’^ (F)*" (3.17) 

is continuous for any /3 G M and C £ C. Bowever, it has better properties under additional 
assumptions described in terms of the operator pencil 

C 9 7^2 SIa (C,?? 2 ) = (x,(Ii -I- <,<92 -7 7772 ) - A,iV° (x,(9i -7 <,^2 -k im)} ■ (3-18) 
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which corresponds to the problem (I2.9p . (12.101) . (12.6p with the fixed parameter A = G M and 
the dual FBG-variable 

r/= (C,r/ 2 ) G Y, (3.19) 

where ( G [0, 27r) is taken from (13.lip . Concerning 7/2 as a spectral parameter, we regard ()3.18p 
as a quadratic pencil in r] 2 , which is a particular case of a holomorphic spectral family, see m 
Ch.l]. 

The following assertion is proved in [33], see also |39l Thms 3.4.7 and 5.1.4]. 


Theorem 2 The operator (3.11) is Fredholm if and only if the segment 

T('/3 = {??2 G C : ReT]2 G [0, 27r), Imr /2 = P} 


(3.20) 


in the complex plane C is free of the spectrum of the pencil If a point of the spectrum 

belongs to T^, then the range of the operator not closed. 


Let us assume that = A*^ satisfies (|2.19l) and, in particular, that the segment Tq is free of 
the spectrum of 2l;s,ti (C)') for all C G [0, 27r) . Indeed, if r ]2 G Tq belongs to the spectrum, then 
At* = A (((^, 7 / 2 ) becomes an eigenvalue of the problem ()2.9p . (I2.10p . ()2.6p and therefore falls into 
some spectral band B^. By the analytic Fredholm alternative, see, e.g., [ 20 I Thm. 1.5.1] or |42l 
Thm.VI.14], we conclude that the spectrum of the pencil (I3.18P consists of a countable set of 
(normal) eigenvalues without finite accumulation points. For 7/2 G C, the spectrum is invariant 
with respect to the shifts 7/2 7/2 ± 27r along the real axis, by the same argument as in ()2.13p 

and (|2.14p . Moreover, it is mirror symmetric with respect to the real axis because, for fixed real 
A** and the problems (|2.9p . (I2.10p . (|2.6p with 7/2 and 772 are formally adjoint. Finally, under 
the assumption (I2.19p . there exists a positive continuous 27r-periodic function 

[0,27r) 9C^/3 o(a»;C) GM+ (3.21) 


such that the rectangle 

|772 G C : Re772 G [0,27r) , llm772l </3o (^A**; c) | (3.22) 

does not contain any eigenvalue of the pencil (C, ■), but the segments 
We further put 

f3o = min|/3o : C G [0,27r)| > 0. (3.23) 

When A** and f are real, the problem (13.lip . (|3.9I) in the infinite strip IT is formally self- 
adjoint, which just means the validity of the Green formula 

( x , di + iC, 82 ) [/ - A»[/, y) ^ (x, di + iC, 82 ) u, g) ^ (3.24) 

= (u,L^ ix,8i+iC,82)V - + (u,N^ (x, + <, ^2) 

for all t/, G G L^perj (n)”- Hence, our assumption ()2.19p and Theorem [2] assure that the operator 
2o (A**; C) in the Sobolev space (H)” is Fredholm of index zero. The next theorem follows 
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from a general result in [33] (see also [391 §3.4 and §5.1]) and it concerns the finite-dimensional 
subspace 

kerT^ ^ ^/ 3 ,pertt (n)" : U'^ satisfies the homogeneous (3.25) 

problem (13.lip . (13.91) with F'^ = 0, G** = O}. 

Theorem 3 Under the condition Ii2.19\) with At* = the subspace ker (A**; C) is independent 
of the weight index (5 € (—/3o (A**; C) , /3o (A**; C)); where /3o (A**; C) is determined in Ii3.21]} - li3.2^) . 

Theorem 4 Let A** satisfy 112.19\) and let (5 € (—/3o (Att;C) )/3o (-^^sC)) where Pq (Att;C) is taken 
from \3.21]) . For any fixed C £ [0,27r), the problem \3.11\) . i3.9i) with the right-hand side 

{f«, G“} G (n)” X (r)” (3.26) 

has a solution [/** G kk|pertt (n)"’? if o,nd only if the compatibility conditions 

= 0, VyGkerr_^(A»;C) (3.27) 

is met. This solution is defined up to an addendum in kerTo (Att;(') = kerT±^ o,nd, if it 

is subject to the orthogonality condition 

(u\v)^ = 0 , vy G kerTo (a«;C), (3.28) 

then it becomes unique and meets the estimate 


uh (n) 


<c(/3)(||T#;T2(n) 


+ 




(r) 


(3.29) 


In the case |/3| < Po{^'^) the constant c(/3) can be chosen independently of f € [0,27r). 


These theorems mean that any solution [/** G Hper^ (n)” of the homogeneous problem 
(13.lip . ()3.9I) has exponential decay at infinity. Moreover, such solutions form the set of all 
defect functionals (I3.27P of the operator T^ with |/3| < fio (Att;C) • Finally, assume that 

G Ffpgrj (n)" and G ^ solutions of the problem 

(I3.1ip . (|3.9p . and the right-hand side {F‘^,G'^} satisfies (I3.26P for both fi = 0 and fi = Then, 
clearly, and may differ by an element of subspace (I3.25P only, and if the orthogonality 
condition (|3.28l) holds true for both, then 


3.3 The spectrum of the model problem in the strip 

Let us consider the spectral problem in II, which is the homogeneous (T** = 0, G** = 0) problem 
(|3.1ip . ()3.9p for the spectral parameter A**. Its variational formulation is: find a number A** G C 
and a non-trivial vector function t/** G (n)” such that 

[a^d (di + <, 82) uKD (di + <, 82) = A# (c/», , vy# G (n)". (3.30) 
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The sesquilinear Hermitian form on the left of (I3.30p is evidently positive and closed as a 
consequence of Korn’s inequality (|1.10p in the finite cells and {x € 11 : \x 2 \ < h}. Thus, the 
problem (j3.30l) is associated [H 10.1] with a positive self-adjoint operator (C) which has the 
differential expression L'^ {x,di + zC, 02 ) and the domain 

V (o) = {t/« G (n)" : iv» (x, + iC, ^2) 1/# (x) = 0 , x g r}. (3.31) 

According to Theorem [2j the essential spectrum of the operator (C) and therefore of problem 
(I3.30p equals 


Ue, ( a ^ (O) = [JbI (C) , bI (C) = {Afc ( 7 ?) : r? = (C, m) e Y} , (3.32) 

fcGN 

while according to (j2.15p 

Bt= U (3.33) 

Ce[ 0 , 27 r) 

If 

A« ^ bI (Co) , Vfc G N, (3.34) 

= dimkerTo ^A^;Co^ > 0, (3.35) 

hold for some Co G [0,27r), then A^ is an eigenvalue in the discrete spectrum a^i {A^ (Co)) and 
kerTo (A**; Co) is the corresponding eigenspace. By the continuity of functions (|2.13l) . we have 

X^BliC), VA^gN, Ag (a#-5o,A« + 5o) , Cg (c“-eo,C^ + eo), (3.36) 


where eo and Jq are positive and the points C and C — 27r are identified due to the evident 
27r-periodicity. Hence, by general results of the perturbation theory of linear operators (cf. [231 
Ch. XIII], [2S1 Ch. 9]) the point (C^ A**) G [0, 27r) x 1R+ is the intersection of k'^ continuous curves 
A = A^ (C), which can be extended either periodically onto the whole semi-interval [0,27r) 3 C 
or have endpoints (cf, ) at the edges of the spectral bands (1021) . see Fig.E Thus, we have 
obtained an at most countable family of bounded, connected and closed sets, that is, segments 
with endpoints included, 

Bl = {X = Xk (C) : C G [Cfe-,Cfc+] C [0,27r]} C M+. (3.37) 

The segment (I3.37P may be covered by a spectral band (I2.15p = (|3.33l) . or may intersect a spectral 
gap and in this case be contained in the additional part (11.251) of the essential spectrum aes (A.) 
of the problem (ll.lSp , ()1.19p ; this fact has been mentioned in Section 11.41 and will be of interest 
later. 

3.4 The regularity field for the open waveguide 

We aim to verify the formula 

aes (Al) = aes (aI°) U u’*, (3.38) 
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Figure 5: Graphs of eigenvalues of the problems in and 11 with continuous and dotted lines, 
respectively. The spectral bands are the boldface vertical segments on the A-axis. 


and to this end we take a point A G M+, which does not belong to set (13.381) . Here, (|3.38p 
contains the essential spectrum (I2.16P of the purely periodic problem (11.41) . (11.51) and also the 
union (cf. formula (11.261) ) 

U (3.39) 

CG[ 0 , 27 r) 

of the discrete spectra of the operator family 27 r) ’ P-3ip . As explained above, the 

set (j3.39p consists of a union of at most countably many segments (j3.37p . 

Let the inhomogeneous problem (11.181) . (11.191) . namely 


L (x, V) u (x) — Xu (x) = / (x), X G n, 
N (x, V) u{x) = g {x) , X G 912, 


(3.40) 


have the right-hand sides 

{/, g} G {nf X 77^/2 (912)" . (3.41) 

We introduce a cut-off function Tq G C°° (M^) such that 

Xq (x) = 1 for either |x 2 | > h + d, or xi < —1 — d, (3.42) 

Xq (x) = 0 for |x 2 | < h and xi > —1, 0 < Tq < 1 . 


The function Xq is equal to 1 outside the semi-strip (—1 — d, -foo) x[—h — d,h + d) but vanishes 
inside the smaller semi-strip (—l,-|-oo) x {—h,h) which contains tu^ and thus also the open 
waveguide. Putting 

f = To/, 5° = Ao5 (3.43) 

gives us vector functions defined in 12^ and 912*^, respectively. Moreover, 


f-L^ (fl°)|| < II/; (fl)||, 


5 °; 2 / 1/2 (^^ 0 ) 


< 


5 ;// 1/2 ( 9 fl) 


(3.44) 


The first inequality (j3.44p is evident, while the second one is a consequence of the following 
observation: in view of (|3.42l) and (11.31) the function Xq equals either one or zero on each 
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connected component of the boundary This also shows that the commutator Tq] will 
be null the in second formula of (j3.49p . 

Since A stays out of the set (I3.38P by assumption, the condition (|2.19p is met and, according 
to Section ESI the problem (I2.19P gets a unique solution G (id®)”' with the estimate 


\\<c (||/°; (0°)II + 5°; {dn^) 


< 


< c 


;L^{n)\\+ 


(3.45) 


This solution forms the first component in the representation (notation will be introduced later 
step by step) 

u = n{X) {/, g} = Xqu^ + (3.46) 

of a parametrix for the boundary-value problem (|3.40p : a parametrix is by definition a contin¬ 
uous operator 

7^ (A) : {VLf X ^ ^2 

such that the mapping 

{L, iV} 7^ (A) - I : (0)^ x {d^f {Vif x {dXlf (3.48) 


is compact, where I stands for the identity. 
We have 


L {Xou^) = Xof + [L, Tfo] = Xif + [L^ Tfo] , 
N {Xou^) = Xog^ = X^g, 


where the term with the commutator [L*^, Tq] admits the estimate 

\\[L,Xo]u^-,H^ (f2)|| = \\[L^,Xo] u^-H^ (f20)|| < 

<c(||/;L2 {n)\\ + \\g-,H^/^ (5f2)||) 


(3.49) 


(3.50) 


The difference = u — is to be sought from the problem 

L{x,V)u^ (x) = {x) = (l - Xoix)"^^ f {x) - [L,Xo]u^ {x), x £ Q, (3.51) 

N (x, V) (x) = g^ (x) = ^1 — T’o (x)^^ g (x), x € dfl, 
where the supports of the right-hand sides are located in the closed semi-infinite strip 


[—1 — d, -t-oo) X [—h — d,h + d] . (3.52) 

We now introduce the cut-off function 


-^tt (^) = x{xi- R), (3.53) 

where x and R are taken from (j2.3np and (|1.13p . respectively. The products 

P = X^f\ g^ = X0^ (3.54) 
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are defined in the domain see (13.31) . and its boundary dQ^. Moreover, the estimate 


/«;i?(!!«) 

+ 


<c(||/^l2 (0)11 + 

(50) 

) (3.55) 




<c(||/;L2(0)|| + 


) 



is valid with any weight index /3 £ R, because < Cj^ in the strip (I3.52p . where the supports 
of and are contained in. 

We now make use of our assumption that A does not belong to the set (|3.38p . As a con¬ 
sequence, the operator Tp{X;() of the problem (I3.11jl . (13.91) is an isomorphism between the 
weighted spaces in (I3.17P for all 


C £ [0,27r) and f3 £ [0, /3o (A)), 


(3.56) 


where /3o (A) > 0 is determined in (j3.23l) . In this way we apply the partial FBG-transform to 
the boundary-value problem (j3.4p with the right-hand sides (|3.54p : this yields a problem of the 
form (|3.11l) . (13.9p . and we find a unique solution U'^ £ (n)” for it. The above mentioned 

isomorphism property of (A; () guarantees the estimate 


p27V 

Jo 


(n) 


/ 2 2 \ 

(^F#;L^(n) + Gi-wy\r) jdC (3.57) 

fhLl 


< c 




where the latter inequality is based on ()3.16l) (and ()3.4p ). We then employ the inverse transform 
(j3.10p and obtain a solution u'^ £ together with the relation 




< Cy3 


J-,wl 


< Cf 


I'2'k 

L 




dC. 


(3.58) 


As a result, we have defined the second component in representation (I3.46p . 


Remark 5 The same cut-off function US. 53]) multiplies in i3.46\ ) as well as and in 
jS.Sfi )- In Section \4^\ we explain how the second term in ^3.46\) must he corrected in the case 
of the L,V- and Y,Z-shaped waveguides (see Fig.\^ and Section lTfR) . 

The estimates (j3.58p . (I3.59p and (I3.55h show that the product in (I3.46h cannot prevent 
the continuity of the parametrix ()3.47p . Furthermore, we have 


L 

N 


«) = + [L, Aft] J = + [l#, a# 


u" 


(3.59) 


2„1 


= Aj5 » = XI g 


where 


[L,Xi]J-,W^ <c <c(j\f;L^n)\\+ g-,H^/^{dn)y (3.60) 
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Since the supports of the coefficients in the operator are located in the vertical strip 

{x : xi e [i2, R + d], X2 G M} , the exponential weight in the norm on the left-hand side of (|3.60p 
helps to prove that the mapping 


L^inrxH^/^ {dnr3{f,g}^ 




(3.61) 


is compact. Indeed, the mapping can be presented as a sum of a compact operator (due to the 
compact embedding C in a bounded domain, the rectangle [R,R + d] x [—t,t]), and a 
small operator with norm O (due to the weight which grows exponentially as |x2| > t 

and t —)■ -l-oo). 

All terms with the compact embedding property, e.g. (j3.6ip . can be excluded from forth¬ 
coming considerations. Hence, owing to the inequality (I3.50p and formulas (I3.58p . (I3.59p . it 
remains to deal with the problem (|3.40l) with compactly supported right-hand sides 

/' = (1 - ^2) (1 - Xf) f, g^={l- Xi) (1 - Xl) g. 

Since this problem is elliptic, recall Section 11.21 classical results in mm give a vector function 
G (^)” such that 

L [X^v^) = X^f^ + [L, Ab] = f + [L, A,] u\ 

N (xm^) = Afe/ = g\ 


(H) 

-b 

[L,Afe] (H) 

<c( 

(H) 

-b 

g^-H^/^ {dn) 

)s 




<c(||/;L2 (H)|| 

-b 

g-H^/^ (90) 

) 



where the cut-off function A), can be chosen as 

( 3 ^)=n n ^ ^ 

± p=l,2 


and p G N is sufficiently large. 

So we have constructed the parametrix TZ (A) with all necessary properties, see (j3.47|) and 
(j3.48p . Since the problem ()3.40p is formally self-adjoint, we also have proved the following 
assertion which, in particular, shows that the given point A belongs to the regularity field of 
the operator A of Section 11.31 


Theorem 6 Assume A G M+ \ Ues (Al), c/. i3.38il . The homogeneous problem Ii.5l) has 

a finite-dimensional space ker {A — A) of solutions in ■ The inhomogeneous problem 

^ 3.409 admits a solution u G (^)”; if and only if the right-hand side i3.41\ ) satisfies the 
compatibility conditions 


+{g,v)g^ = 0, VuGker(^-A). 

This solution is defined up to an addendum in ker (Al — A), and if in addition the orthogonality 
conditions 

{u, v)q = 0, Vu G ker (Al — A) 
hold true, then it becomes unique and meets the estimate 


u;H2(0)|| <c(||/;L2(0)||+ g-H^/\dn) 
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3.5 The essential spectrum of the open waveguide 

We now conclude with the main result of the paper. 


Theorem 7 The essential spectrum of the operator A of the problem /11.19\) equals 

aes{A) = aes{A^) U see US. 38\) . where the components aes {A^) and are defined in 
Seetions \2.S\ and \3.3\ bv the formulas i2.16\) . I13.39\) . respectively. 


Proof. Thanks to Theorem [6l it suffices to construct singular Weyl sequences for the 
operator A at all points in the set (13.381) . If A G cJes , we may take such a sequence from 
Section [231 because the support [2-^ , x [2-^ , 2'^+^] of the cut-off function (I2.31|) and the 
entries (I2.32p of the sequence do not touch the semi-strip {x : xi > 0, |x2| < h}, where the open 
waveguide lies in. 

Let A belone: to the interior of segment (13.371). and recall that the endooints live in a es (AI^) , 
if they exist. By definition, there exists a C S [0,27r) such that the problem (13.lib . (13.9h admits 
a non-trivial solution U'^ G -^pertt (n)”i generates a Floquet wave in the xi-direction 

u^x) = e<^^U^ {x) (3.62) 


satisfying the homogeneous (/^ = 0, <7^ = 0) problem ()3.4p in the periodic domain ()3.3p . 
By the constructions in Sections 11.31 and 13.11 this domain coincides with Tl in the half-plane 
{x : xi >0}, where the matrices A and become equal to each other, see (|1.16l) . (I2.33p . (13.11) 
and (11.151) . (13.21) . We localize the wave (I3.62h by a cut-off function similar to (I2.31h . namely 


{x) = Xj (a^i) X {x2 + 2^+^) y (2^+^ - X2) 


(3.63) 


where Xj and x are taken from (|2.29p and (j2.30p . The function Xj vanishes outside the rectangle 
[2-^, 2-^+^] X [—2-^+^, 2-^+^] and we will deal with indices j G N such that 2^ > R, cf. (I1.13h . We 
set 

(x) = X - (x) (x), (x) = (n) (x) (3.64) 


and observe that, according to our choice of cut-off functions, both vector functions in (I3.64h 
satisfy the boundary conditions (|1.5p . 

Recall that by Theorem [3] the function u^ (x) decays exponentially as X2 ±00, so we can 
compute the (n)-norm of as follows, cf. (I2.33p . 


(52) 


2J+1-1 1 „2J+l-l 


> 


E 

j=2i+l 



0 X_2J+i+l 


uH 


X 


dxidx2 > 


(3.65) 


> (2^+^ - 2 ^ - 2) (n) — C exp (-2/32^+^)^ > C{t2^ ctt > 0. 


We also obtain 


fij ^ lyij _ XyiJ = - XyiJ = \L\ 


u" 


and notice that the supports of the coefficients in the commutator 


lKxI 


(3.66) 

belong to the union 


of two horizontal and two vertical rectangles of sizes d x 2 ^ and x d respectively (they are 
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overshadowed in Fig.^a). Due to the exponential decay of the horizontal rectangles only 
cause an infinitesimal input into the (D)-norm of (I3.66P as j —^ +oo, and the input of the 
vertical ones stays uniformly bounded in j. These mean that 




(D) 


/«^ ;L2 (D) 




An application of the Weyl criterion finishes the proof. ■ 


Remark 8 As was commented in Section \1.4\ set 113.39\) in the representation i3.3^) of 
(Tes (-4) can have a non-empty intersection with the essential spectrum a^s of the problem 

in the purely periodic perforated plane 


4 Possible generalizations of the results 

4.1 Concrete problems in mathematical physics 

(i) Scalar equations. Let D (V) = V and let A be a Hermitian, positive definite 2 x 2-matrix 
function possessing the properties described in Section 11.31 Then n = m = 2 and 

L{x,V,) = -V^A{x)V (4.1) 

becomes an elliptic second-order differential operator in the divergence form. The algebraic 
completeness is evident, and in (jl.9p we have pv = 1- We consider the following boundary 
conditions and problems: the Neumann condition 

n {x)~^ A{x)'Vu{x) = 0, X € di}, (4-2) 

which is nothing but (|1.5p with the co-normal derivative (II.8p . and the Dirichlet boundary 
condition 

u{x) = 0, X G dfl, (4.3) 

for the generalized Helmholtz equation 

— V~’'A (x)Vu (x) = Ap (x) u (x), X G fl, (4.4) 


and the Steklov spectral problem 


—Au(x) = 0, X G Q, (4.5) 

di,u (x) = Xu (x), X G 9D, 


where the spectral parameter A appears in the boundary condition while = v (x)"*" V and 
A = V^V denote the outward normal derivative and the Laplace operator. The problem (14.41) . 
()4.2p occurs in acoustics, and the problem (j4.4jl . (14.3p with L = A and A (x) = I in the theory 
of quantum waveguides, see Remark 0 and [71 [HI [101 [121 [T3] . Moreover, (|4.5I) is related to the 
linear theory of water waves, cf. [9l[30] and Remark 1101 
(ii) Elasticity. Let n = 2, m = 3 and 


D(V) 


di 0 2-1/292 \ 

0 82 2-1/29, J ’ 


A. 

dXn 


1 = 1 , 2 . 


(4.6) 
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Interpreting u = {ui,U 2 )~^ as a displacement vector, we employ the Voigt-Mandel matrix nota¬ 
tion in elasticity and introduce the strain and stress columns 

e{u-,x) = (^eu{u-,x) ,e 22 {u;x) , 2 ^/‘^ei 2 {u;x)^ =D{V)u{x), (4.7) 

<T {u; x) = A (x) £ {u; x) = A [x) D (V) u (x), 

which are composed from the Cartesian components £jk (u) = 2~^ (djUk + d^uj) and cTj^ {u) of 
the strain and stress tensors, respectively. Here, A [x) stands for the Hooke matrix of elastic 
moduli, and it is real, symmetric, uniformly bounded and positive definite. The matrix D (V) 
in (I4.6jl is algebrically complete and qd = 2 , 'm (II.9p . see [iQl §3.7.5] and, e.g., [Ml Example 
1.12], m- We assume that A (x) meets all requirements listed in Section 11.31 

Time harmonic elastic waves with frequency k > 0 satisfy the system of differential equations 

— diCTji {u\ x) — 820 -j 2 [u] x) = (x) Uj {x ), j = 1, 2, x G H. (4.8) 

If p (x) = const > 0 and X = system (14.8p . in view of (|4.7p and (j4.6|) . takes form (II.4p while 
the traction-free boundary condition on dQ reads as (|1.5|) . 

The problem on the oscillations of a homogeneous (constant A), perforated elastic plane is 
surely interesting for the engineering applications. To include composite elastic materials into 
our considerations, we must deal with the variable piecewise smooth Hooke matrix A (x) and 
material density p {x). Notice that the same modification can be applied to (j4.4|) as well. In 
the next section we will show how to get rid of the smoothness assumptions made until now 
when adapting our method to the variational formulation of the elasticity problem. Finally, as 
a possible application we also mention the Dirichlet problem (14.3p for the elastic displacement 
vector ti, which has in the two-dimensional case a clear mechanical interpretation meaning that 
the boundaries of the holes in a thin elastic plate are rigidly clamped. 

{in) Piezoelectricity. We set n = 3, m = 5 and 


Z)(V)^ 


©2X2 \ p 

Oix3 D^(V)^J’ 


1 0 0 
0 1 0 
0 0 0 


diag{l, 1,0} , 


(4.9) 


where the 3 x 2-block (V) is taken from (|4.6p and (V) = V. The column u = 

consists of the displacement vector = {ui,U 2 )~^ and the electric potential = U 3 so that the 
superscripts M and E indicate mechanical and electric fields, respectively. The so-called smart 
piezo-devises are able to couple these fields of different physical nature, and this phenomenon 
is described by the following system of three differential equations, see [laiiiiES], 



D{-V)'^ A{x)D{V)u{x) = Xp{x)Eu{x), x £ Q. (4.10) 


Here, p {x) > 0 is the material density and the matrix A {x) is written blockwise as 


A{x) 


(x) —A^^ {x) ) 


(4.11) 


where [x] denotes the elastic Hooke matrix, A^^ {x) the dielectric matrix and (x) = 
^EM 

(x)"*" the piezoelectric moduli. All matrices are real and (x) and A^^ (x) are symmetric 
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positive definite, hence also the matrix (14.lip is symmetric, but it is not positive definite due 
to the minus sign in the bottom right-hand block. This reflects the intrinsic transformation of 
the elastic energy into the electric one and vice versa in a piezoelectric body. At the same time, 
the electric potential does not affect the kinetic energy at low and middle frequencies and 
therefore it is absent on the right-hand side of (14.101) . cf. structure of the diagonal matrix E in 
(j4.9p . We emphasize that in spite of the minus sign in (14.lip the piezoelectricity system ()4.10p 
is elliptic in the Douglis-Nirenberg sense (cf. [351 Example 1.13]). 

In Section 14.31 we demonstrate a reduction scheme from [36] , which allows us to apply the 
above results to the piezoelectricity system with various boundary conditions. 

{iv) Plates. Let n = 3, m = Q and 


Z1(V)^ = 


0 2 - 1/292 

d 2 2-i/29i 
0 0 


0 0 0 

0 0 0 

d1 d1 21/2 02 



The Dirichlet problem (j4.10p . (14.3p with matrices (j4.12p and 


©2 


x3 


V ' (V) 


(4.12) 


A{x) = 


24BM ^BB 


E = 



(4.13) 


describes the Kirchhoff model of an anisotropic inhomogeneous plate with rigidly clamped 
boundaries of holes, which means 

rt (x) = 0 e d^us (x) = 0, x € dQ. (4-14) 

The vector u = (tti, rt 2 , Ms)"*" includes the longitudinal displacement (ui, U 2 )~^ and the deflection 
U 3 . The matrix A (x) is real symmetric and positive dehnite, but in the case of elastic symmetry 
it becomes block-diagonal, i.e. = (A^''^)"'" = Osxs- In this case the Douglis-Nirenberg 

system (|4.10p decouples into two second-order equations and a fourth-order equation. In Section 
I4.3l we show that our scheme works also in the case of high-order differential equations including 
the Kirchhoff plate (see also m)- 


4.2 Operator formulation of the variational problem 

The integral identity ()1.2np corresponding to the spectral problem (jl.lSp . (|1.19p makes sense 
even in the case the matrix A and scalar p are just bounded, measurable, and the boundary 
9D is Lipschitz. It is important that all results on the solvability of model boundary value 
problems in Section 14.11 are easily adapted to their weak formulations in the Sobolev and Kon¬ 
dratiev spaces (see m Ch.2], and for the periodic case [37|). It is straightforward to include 
into our consideration any type of boundary conditions, which are covered by the symmetric 
Green formula (cf. [31[ §2,2]), for example, the Dirichlet conditions (|4.3p or mixed boundary 
conditions. 

Remark 9 Lipschitz domains occur for example in the grating of quantum waveguides with 
long-haul thickening as in Fig.\^ a). Notice that the domain of the operator A differs from 
(D) n Hq (D) , cf. and, e.g., ffPl Ch.2]. 
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Figure 6: Other periodic geometries suitable for the present method. 

If the density p is not a constant, it is useful to change the operator formulation of the 
problem 

{AD{V)u,D{V)v)^ = X{pu,v)^, (4.15) 

in order to apply the theory of self-adjoint operators in Hilbert space. Namely, having in mind 
the Korn inequality (11.171) . we introduce in the specific scalar product 

{u,v)^^ = a{u,v;n) + {pu,v)f;^, (4.16) 

and then dehne the continuous, positive, symmetric, hence self-adjoint, operator 5 in by 

{Su,v)^^ = {pu,v)^, Vu,uGio. (4.17) 

This turns the problem (|4.15p into the abstract equation 

Su = in Sj (4.18) 

with the new spectral parameter 

? = (1 + A)-^ (4.19) 

The above-mentioned theory readily applies to the equation (|4.18l) . 

Remark 10 In the perforated plane 11 (where is not empty in and also the 

Steklov problem |.^.5[ ) reduces to the equation with parameter i4-19\ ). According to 

we endow in this case the space (H) with the scalar product {u,v)^ = {Vu,Vv)q + {u,v)qq 
and define the ’’trace operator” S by {Su,v)^ = {u,v)qq . 

4.3 Reduction to integro-differential equations 

Concerning the piezoelectricity problem with matrices (j4.11l) . (14.9p or the plate problem with 
matrices (j4.1.Sp . (I4.12|) . we hrst mention that the results of the papers [33l [M] can be applied 
here, since they deal with general boundary value problems for Douglis-Nirenberg elliptic sys¬ 
tems. However, the degenerate matrix E on the right-hand side of (Id.lOj) hampers the use of 
the theory of self-adjoint operators in Hilbert space. 

For the perforated Kirchhoff plate {il^ 0), the above mentioned trick works with the new 
scalar product and operator 5, 

(u, v)^ = {AD (V) u,D{V) + (pEn, Eu)^^, (4.20) 

{Su,v)^^ = {pEu,'Ev)^ , yu,v&S) 
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where matrices are taken from (j4.13l) and (I4.12p . Indeed, owing to the Dirichlet clamping 
condition (I4.14p the bilinear form (I4.20p is a scalar product in = Hq x Hq [dO) . 

The piezoelectricity problem (j4.10l) . (II.5p requires a much more elaborate process. Following 
[36], we reduce the corresponding variational problem 


^ (V) U^, (V) 

+ (V) u^,D^ (V) - (a^^D^ (V) u^, (V) 

to 

where 

{V)Ru^,D'^ 

and = Ru^ G R^ (fl) is a solution of the scalar Neumann problem 

Here, the space R^ (Q) is defined as a completion of (H) in the norm 


Vu^;L^(n) 


+ 


u^;L^ (K) 


1/2 


(4.21) 

(4.22) 


(4.23) 


and K is a compactum in H of positive area. Since constants fall into R^ (Q ), the problem 
p4.23p includes one compatibility condition, which is obviously met because the functional on 
the right-hand side of (I4.23|l degenerates on constants. In this way, the bilinear form (I4.22h 
is well-defined, symmetric and positive, see [36| for details. Now the trick with a new scalar 
product in (fl)^ applies again. 

Similar modifications work for other types of boundary conditions in piezoelectricity and 
plates as well. 


4.4 Other geometries 

To simplify the notation we have always used the covering of the plane with unit squares. The 
cells can as well be rectangles or parallelograms, because an affine transform does not change 
the crucial properties of the matrices D (V), A{x) and operators (II.7p . (II.Sp . We mention 
that according to [29|, the elasticity problem in Section 14.II (ii) preserves its matrix form under 
any affine change of coordinates, though some special (non-physical!) columns of strains and 
stresses will then appear. A similar procedure applies to piezoelectricity and plates in Section 
lO {iii) and (iv). 

Other types of planar coverings can be treated in a similar manner, for instance, the diamond 
and honeycomb shapes, see [l5|. We do not touch upon this generalization, which would require 
a total modification of the notation. 
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Figure 7: a) Perforation of another type, b) E-shaped open waveguide. 


The purpose of the restriction (11.31) was solely to simplify the definitions of cut- off functions 
in (j2.29p - (|2.32j) etc., but other settings of holes, cf. Fig.[3a), can be treated by our method as 
well. We remark here that the homogeneous boundary conditions (|1.5p and (I1.19P are naturally 
included into the definitions of the domains (jl.lip and (jl.2ip of the operators and A, and 
hence one has to take care that multiplication with the plateau functions (I2.3ip and (13.631) does 
not spoil these conditions, otherwise the arising discrepancies must be compensated. A simplest 
way to avoid the discrepancies is to construct cut-off functions, which respect the geometry and 
surround all holes, like indicated in Fig.Qa), where the support of |VAj| is shaded: no serious 
change of calculations would follow. Of course, no modification of cut-off functions is needed 
in the case of an intact periodic medium and Dirichlet conditions (14.3p . 

Our methods are sufficiently general to study the essential spectrum of the problem (jl.lSp . 
(|1.19p in the layer {x = {y,z) : y = (^ 1 ,^ 2 ) € 1^1 <77} perforated periodically and per¬ 
turbed inside the semi-infinite cylinder M_|_ x {—h,h) x see Fig.[6lb). As in two- 

dimensional case, in this three-dimensional problem there again appear only two model spectral 
problems, one in a cell and another one in an infinite perforated prism, which can be exam¬ 
ined with results in [451 [2H] and |33l EH] . In particular our approach applies to elastic infinite 
three-dimensional plates, which are perforated or have periodic bases. On the other hand, 
our approach does not yet help to examine a problem in the space with a triple-periodic 
perforation and with an open periodic waveguide inside a semi-infinity circular cylinder. 

4.5 Joints of open waveguides 

The waveguides depicted in Fig. [3] must keep the periodicity along all their branches so that 
axes of inclined branches in Fig.|3l a), c), cross the xi-axis at angle (p with a rational tan cp. Then 
each branch gives rise to its own model problem of type (13.lip . (13.91) in a perforated strip Ilg 
which is perpendicular to the branch axis and may have width different from the main period 
1; here q = 1, ...,Q and Q is the number of branches. In this way the essential spectrum of 
problem (jl.lSp . (|1.19l) for the joint of open waveguides becomes the union of aes and the 
sets a\, defined in (I3.39P through the discrete spectra of the model problem in Hi, ...,11^ 

respectively. Theorem [7] remains true with this modification, but let us next comment the minor 
changes required for its proof. 

First, let us reconsider the Weyl sequence in Section [231 Since there exists an unbounded 
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angular opening between adjacent open waveguides, we can choose for any j € N a. point 
such that the support of the plateau function 


(x) = Xj -Pi) Xj [x2 - Pi) 


(4.24) 


intersects neither the open waveguides, nor supports of Xi, Using (j4.24l) . the entries 

(j2.32l) of the Weyl singular will have all the properties listed in Section [231 Similarly, in Section 
EE we shift the ’’center” of the plateau function (I3.63|l along the branch axis and so obtain the 
Weyl sequence elements (I3.64p . 

Second, the construction (j3.46p of the parametrix (j3.47p now involves the solutions of 

(|3.4p in the domains with the periodicity ’’cell” 11^; here q = 1, These solutions should 
be located near their branches with the help of some cut-off functions However, the def¬ 
inition (I3.53P does not work properly, since suppT’jJ^ may intersect other branches of the joint 
and therefore new discrepancies may appear. To avoid these discrepancies one may place the 
plateau function between two neighboring branches as indicated in Fig. [3 b) by overshading. 


Since supp 




is located in a p-neighborhood of the sides of an angle domain, the expo¬ 


nential weight in the Kondratiev space still leads to compact mappings of type (I3.6ip and as 
a consequence the mapping ()3.48p remains compact. Other steps in our proof of Theorem [3 
remain unchanged. 

We finally mention that in our notation the E-shaped joint of open waveguide of Fig.jTlb), 
only has one branch, and therefore it is directly covered by Theorem [3 
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